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sions and numerical diagonalization. We are mainly interested in results for general chiral 
angles but we consider also the superintegrable case. For the parameter values consid- 
ered, we find that the band structure of the low-lying part of the excitation spectrum 
has the form expected from a quasiparticle picture with two fundamental particles. 

Studying the iV-dependence of the spectrum, we confirm the stability of the second 
fundamental particle in a limited range of the momentum, even when its energy becomes 
so high that it lies very high up among the multiparticle scattering states. This is not a 
phenomenon restricted to the superintegrable line. 

Calculating a non-translationally invariant correlation function, we give evidence that 
it is oscillating. Within our numerical accuracy we find a relation between the oscillation 
length and the dip position of the momentum dispersion of the lightest particle which 
seems to be quite independent of the chiral angles. 
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1 The TZq— chiral Potts model 



Recently, the Z n -symmetrical chiral Potts model [p], [2|] has attracted particular interest 
Pf| — fllT] [] because of its special integrability properties and the rich structure of its phase 
diagram. In this paper we consider the ^-quantum chain version of the model |L2] . 



The Hamiltonian of the model contains complex coefficients and this gives rise to a large 
number of ground state- (and also higher level-) crossings, and to the appearance of 
both massive and massless oscillating phases. We present evidence that the spectrum 
has a quasiparticle structure all over the massive high-temperature phase. We also study 
the peculiar momentum dependent effects which appear due to the parity-noninvariance 
of the model: A quasiparticle may overlap with a two-particle spectrum for a limited 
range of its momenta. We investigate whether by immersing into the continuum the 
particle remains stable or ceases to exist as a quasiparticle. Such an effect has been 
discovered recently in || for the superintegrable case of the chiral Potts chain. Finally, 
we present finite-size calculations which establish directly an oscillatory behaviour of 
certain correlation functions in the massive high-temperature phase. 

The Z 3 -chiral Potts quantum chain with N sites is defined by the Hamiltonian: 



o N 

(3) 



-== £ {e~ i ^a j + e^/Vt + A (e^r^ + e^ 3 r+r j+1 ) } . (1) 
V3 i=1 



The operators <jj and Tj act in a vector space d 3 located at site j. They satisfy the 
relations 

ct.V, l> ; ^/ : a J = T] = l; where u = e 2 ™ /3 . (2) 



We shall use the representation 
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In this paper we will consider the ferromagnetic case < A < oo. Using periodic 
boundary conditions Tjv+i = Ti, the Hamiltonian ([!]) commutes with the translation 
operator T and with the Z 3 -charge operator Q = YljLiCj- Their eigenvalues have the 
form e lP and respectively, leading to conservation of the momentum P (which for 
finite chains takes the N discrete values < 2irk/N < 2tt with k = 0, 1, . . . , N — 1), and 
to the conservation mod 3 of the charge Q = 0, 1, 2. See e.g. [p], for the definition of 
the chiral Potts quantum chain for general Z n which for n = 2 reduces to the standard 
Ising quantum chain. 

The Z 3 -model contains three parameters: A which plays the role of the inverse 
temperature, and <fi, <p which are two chiral angles. For (ft = = |, the Hamiltonian ([!]) 
is 'superintegrable', i.e. it fulfills the Onsager algebra |], ||. For 

cos ip = A cos (4) 

(HD has been shown to be integrable due to a new type of Yang-Baxter equation which 
contains two rapidities not enjoying the usual difference property. For generic <f) = <p the 

1 For a more complete list of references see e.g. 
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Hamiltonian (|l|) is self-dual but probably not integrable. The self-dual model is parity- 
invariant only for = <p = 0. In this case ([!]) is just the standard WM^-symmetrical 
Z 3 -Potts quantum chain. The phase diagram of (P contains four distinct phases, two of 
which are massless incommensurate phases. For recent reviews summarizing the known 
properties of the model, see ||, |9|. 

In this paper we shall study properties of the massive high-temperature phase of the 
model which occurs for small A: < A < A c (0, </?). The boundary of this phase, A c (0, <p), 
is known exactly only for a few special values of <f), (p: For <p — (p, it is known that 
A c (0, 0) = 1 (standard Potts case), A c (|, f ) = 0.901292 . . . (superintegrable case) [JJ and 
for (j> — ► 7r , ip —>■ ir we have A c — > ||12|| . 

In Sec. 2 we start explaining our tools for analyzing the spectrum of ([I]), summarizing 
known formulae and presenting simple new perturbation formulae. In Sec. 3 we show 
that the finite-size data follow the band structure expected for a quasiparticle spectrum. 
In Sec. 4 we study the immersion of the Q = 2-quasiparticle into the many-particle 
continuum for a limited range of P. In Sec. 5 we study a non-translationally invari- 
ant correlation function and establish its oscillatory nature which previously had been 
inferred only from wave-vector scaling arguments. A simple relation between the dip 
position of the Q = 1-particle dispersion relation and the oscillation length is found 
to be approximately valid over a large parameter range. Finally, Sec. 6 presents our 
Conclusions. 



2 Methods of calculation 

The complete analytic calculation of the spectrum of (|1|) in the superintegrable case 
cf) = (p = | by the Stony-Brook group (|| and references therein) has established that 
for these parameter values the spectrum has a quasiparticle structure, i.e. in it has 
been shown that in the large chain limit all excitations AEg yr (P, A) above the ground 
state are composed of two fundamental excitations E^P) and E 2 (P) according to the 
following rules: 

m r m r m r 

A%(P,A) = ^E QW (P«), P = ^P« mod 2tt, Q = £ Q« mod 3. (5) 

k=l k=l k=l 

Additionally, the fundamental quasiparticles satisfy a Pauli principle, i.e. 

Q w = Qtf) implies P« ^ P®. (6) 



For more general 0, ip and general lL n numerical results were presented in |10| which 
suggested that such a quasiparticle structure (with n — 1 quasiparticles in theZ„-case) 
is present all over the massive high-temperature phase. Here we want to give further 
evidence for such a structure for the Z3 model (0) for quite general <ft, <p. For this we 
compare multi-particle levels to the basic single-particle levels and check whether the 
former can be composed out of of the latter ones. Single-particle levels are expected to 
show exponential convergence in N to the thermodynamic limit, while composite states 
should show power convergence with the power indicating the number of constituents 

BE3- 



In order to analyze the spectrum, three different approaches will be used: 



1) In each charge-momentum sector we calculate numerically the lowest > 10 energy 
levels of ([!]). For chains of length N < 8 complete numerical diagonalization is easily 
possible, while for up to N < 12 sites we use Lanczos techniques. This way, both the basic 
quasiparticle dispersion relations Ei(P) and E 2 (P) and many continuum levels can be 
obtained. The limit N — > oo may be calculated using standard finite-size extrapolation 
techniques. In general, the convergence in N depends on the mass scale involved. For 
small A the mass scale is large but also for A close to A = 1, where the mass scale is 
small, good information can be obtained as will be shown below. 

2) Compact analytic, though approximate formulae for each lowest level in the Q ^ 
sectors can be calculated for general and — | < if < | using perturbation theory in A. 
We find up to order A 3 : 

E X (P) := AE lfi (P,<f>,<p) 

= 4sin(^^)-i|{^(P)+«-2/3-(« 2 + 2/3 2 )cos0} + O(A 4 ), (7) 

with 

F{P) = (1 - A) cosp + a cos (p + 0) + (3 cos (2p) + B cos (2p + 0) + 2f3 2 cos (3p), (8) 
where A A 

P ~ P ~ "o ' OL — — i-K-w \ i ft ~ 7, 7v\ } (^) 

3 6 cos (—3^) 6cos(|) 

and A = 3/3 2 — 2a/3 + 2a 2 ; B = 2f3 2 + 2a/3 — a 2 . In the superintegrable case we have 
a = (3 = A/(3 v / 3) and A = B = A 2 /9. 

For < if < I the lowest Q = 2 excitation is given by: 

E 2 (P) := AP 2i0 (P, 0, if) = AP 1)0 (P, -0, -if). (10) 

These formulae are expected to be valid for A not too large, i.e. if the distance of the 
level from other multiparticle levels is not too small. Of course, this nondegenerate 
perturbation calculation is useless once a level penetrates into the continuum of multi- 
particle states. We shall see that there are parameter values for which this happens in 
the Q = 2-sector. 

The P-dependence of Ei is contained in T{P) only. So, e.g. to find the value of the 
momentum P = P m i n , for which the lowest energy is reached (because of parity violation 
this is not at P = 0), we have to look for the solution of dJ-"/dp = 0. For non-zero 
0, cos and cos |, we see that p m j n = P mm — | is of order A. So, expanding we get 

Pmin = ^-asin0 + O(A 2 ). (11) 

Observe that, in order to obtain P min to order A, we have to use the second order 
expansion of Ex(P). Starting for A <^ 1 at P mm = |, the position of the minimum 
shifts to smaller values of P for increasing A. Using (|TD|), one obtains immediately the 
minimum of the curve E 2 (P). 

3) For the superintegrable case = ip = | analytic formulae are available from |§ 
for the whole spectrum. We shall use here only the formulae for the energies of the two 



fundamental excitations which are given by 



E l( P r ) = e r = 2 | 1 - A | - / dt t M — - (1 - A)- 
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E 2 (P 2s ) = e 2s = 4 | 1 - A | +- I dt — ,7 ' - (1 - A) 2 , (12) 



^ a ; v 2s (4 vjt 2 - v 2s t + I) fj\ 
vlt 3 + 1 V t 3 - 1 

where the auxiliary variables v r and f 2s are related to the momenta by: 



v(P) 



1 — cos P 



\ 1 - COS (P + 



2?p 
3 ' 



(13) 



/oa /o oa /oa / -^( p ) for < P < 4tt/3; 

v 2s (P) = -v{2n - P) , v r (P) = | + ^ p j for 47r / 3 < 7, ^ (14) 

The integral appearing in E 2 (P 2s ) has a singularity at P = 27r/3, so that the Q = 2- 
quasiparticle level is defined only for 27r/3 < P < 27r. Except for this singularity, there 
is no problem to evaluate the integrals in (0) numerically: using £ = \/£ — 1 instead of 
£ as our integration variable, the integrands are smooth at £ = or t = 1. 

Comparing these different methods, we can keep control about the precision of the 
values to be used for Ei(P) and E 2 (P) for all <p,(p and A < 1. We shall comment on 
such comparisons as we go on presenting our results. Methods 1) and 2) can be applied 
equally well also in the non-integrable and in the non-hermitian situation (the Lanczos 
method can also be adapted to the case of non-hermitian matrices |[L4|| ) which arises for 
complex or ip e.g. if @ is imposed. This way we can cover a large range of 4> and (p. 

As a first example, from which one can judge the usefulness of the different methods, 
in Fig. 1 we show the lowest part of the Q = 1-spectrum for the superintegrable case 
= if = | at the large value of A = 0.90 which is close to the phase transition point. 
We see that the finite-size data for iV = 12 lie still within < 6% on the exact curve 
(|T2{). For small P, even the points for small N move nicely along the exact curve. Of 
course, with only N < 12 we cannot get directly points in the interval < P < 7r/6. 
Considering the large value of A and the vicinity of the phase transition, it is remarkable 
that the 3r<i-order perturbation curve remains close to the exact curve, especially for 
P > 1.27T. It mainly fails to describe details of the dip at small P, and so is not useful 
for the determination of the transition to the incommensurate (IC)-phase (where the 
ground state looses translational invariance). 



3 Quasiparticle- and band-structure of the spectrum 

In order to investigate the band structure expected from the composition of the two 
quasiparticles, in Fig. 2 we start considering the simple situation which occurs for a 
small value of A, here in particular A = 0.1, and a self-dual choice of the chiral angles 
= if = ~ well below the superintegrable line. The single-particle dispersion curves 
can be determined with high precision: the iV = 12 finite-size value for Ei(P) agrees 
to better than 10~ 4 with the value given by formula (0). Also the agreement in E 2 (P) 
between N = 12 and perturbation theory is better than 10 -3 over the whole range 



i 



of P. Both curves are rather flat, corresponding to a large mass in a Klein-Gordon- 
type approximation to the dispersion curve. Of course, at A = the P-dependence is 
completely absent, as is immediately seen from the perturbation expansion. 

From these single-particle dispersion curves we can calculate the shape of the two- 
particle bands in all three charge Q-sectors. These bands are bounded by Eq(P), given 
by 

E+(P) = max {E Ql (P/2 + g) + E Q2 (P/2-g)}, Q = Q 1 + Q 2 mod 3 (15) 

^ — 7r<g<7r 

(and similarly Eq(P) with max 5 replaced by min g ). Now, as long as the P-dependence 
of the quasiparticle dispersion curves has the form 

E Qi (P) = A, + £ By cos (P - aij ), (16) 

3 

the boundaries Eq(P) for two equal particles are obtained for q = 0, tt, respectively, and 
we have simply 

E±(P) = E Ql (P/2) + E Q2 (P/2). (17) 

Of course, we have to use the doubled Brillouin zone for P. In the perturbation formula 
(0) the first term to spoil the form flTfiP is the term ~ A 2 cos (2P — ^) . For the parameters 
chosen in Fig. 2 this is quite small. For two charge conjugate particles (in our Z3-case 
this occurs for Q = 0), with a±j = — «2j and A\ = A 2 , B±j = B 2 j (which is true in order 
A of©,©), we get 



E x (P/2 + q) + E 2 (P/2 - q ) = 2\A 1 + J2 Bij cos (P/2) cos (q - a y ) \ . (1 



so that the band width becomes zero at P = tt. The A 2 -contributions no longer obey 
fully flT5|) with the effect that the lower corner of the zero-width point is rounded off. 

In Fig. 2 we have plotted the nine lowest N = 12-sites levels of each charge and each 
momentum state. As expected in the quasiparticle picture, we see that at all available 
values of P ^ tt there are three bunches of levels with the correct charge filling up and 
following the respective two-particles band ranges. The P-values, where the respective 
bands have zero width is nicely seen to be at P = tt ± 2P min « tt ± tt/6 for the Q = 2- 
and Q = 1-bands, respectively, and at P = tt for the Q = 0-band. For Q — 1 and Q = 2, 
according to whether k is odd or even, the uppermost three, respective two highest levels 
start forming the next bands which at these energies are expected from three-particles 
states. 

Only for P = tt four Q = 0— points appear outside of the bands. The two points 
slightly above E = 8 are within the expected three-particle band but the two points above 
and below E = 6.5 must be two-particle states which, because of power convergence at 
N = 12, are still out of the here very narrow band. Such finite-size effects have been 
discussed in detail for P = in [|llj. There we found that, taking finite-size effects into 



account, for P = the quasiparticle picture explains the data very well. 

In Fig. 3 we show the analogous picture for a somewhat larger value of the inverse 
temperature, A = 0.25, and a (again self-dual) choice of the chiral angles = y? = 67.5°, 
closer to the superintegrable value = <p = 90°. For the Q = 1-ground state the 
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third order perturbation formula still agrees within « 10~ 3 over the whole range of the 
momentum with the N = 12 sites data, however, there are somewhat larger (« 5 • 1CT 2 ) 
discrepancies for the lowest Q = 2-level. 

The band structure in all three charge sectors is of the same shape as before and all 
higher levels lie within the bands (this time we show all higher levels for N = 4, ... 12 
sites but draw only the boundary curves for the Q = 2- band). The approximation ( |P7D 
for the lower bound of the Q = 2-two particle band is no longer perfect, it is interesting to 
see that the finite-size data at the most narrow point P ~ \.22i\ do follow the curvature 
expected. With increasing N the bands are filled densely from the interior which is 
consistent with but of course no proof of the presence of the Pauli principle effect ([]). 

More important is to observe in this Figure that, due to the parity non-invariance, 
the maximum of the E 2 -curve is shifted to the left of P = tt, while the two-Q = 1 -particle 
band is shifted to the right of P = tt. This causes the curve for E 2 to penetrate into the 
two-particle band. This may make the Q = 2-quasiparticle unstable in a limited range 
of the momentum around P ~ 2tt/3. Indeed, the convergence with N of the lowest level 
in the Q = 2-sector deteriorates considerably at these momenta. 

In |9|, |10| it has been pointed out that the Q = 2-particle enters the continuum 
composed out of two Q = 1-particles at P = on the superintegrable line . However, we 
see that for P « 2tt/3 this effect appears already for lower values of 0. With increasing 
A this "threshold" moves to smaller values of 0. It is tempting to suspect that the 
limitation of the range of the momentum for the Q = 2-quasiparticle discovered in || 
for the superintegrable case, see eq. (|T2"D, might be due to this penetration effect. In the 
following Section we try to shed some light on this problem by studying the chain-size- 
dependence of the spectrum. 

Not surprising in view of the penetration of the Q = 2-quasiparticle into the con- 
tinuum, the perturbation formula for E 2 (P) of eqs. (^), fllPl), which for A = 0.25 and 
= ip = 67.5° was still good, becomes quite useless for larger values of A and 0: E.g. 
for A = 0.5 and = tp = 60° it is off by 20 — 25% around P = 0, while in the back- 
ward direction it describes the finite-size data still within < 3%. Going at A = 0.5 to 
= (p = 75°, in the forward direction the A 3 -perturbation formula value is already too 
large by about a factor of three (only around 1.27T < P < 1.5tt the curve comes back to 
the finite-size data within a few percent). 

Looking again into Fig. 1 we see that there are ten points which appear well below the 
curve Ei(P). This happens because close to A = 1 also the Q = 1— quasiparticle enters 
a multiparticle band. This time the band is formed by four of the same Q = 1-particles 
as can be seen by estimating its lower boundary from Ei(P — 3P m i n ) + 3Ei(P min ). So it 
is not trivial that the perturbation formula for Ei(P) remains useful around A ~ 0.9 as 
seen in Fig. 1. 

4 Stability of the Q = 2-quasiparticle in the two- 
particle continuum 

4.1 The superintegrable case 

As we have seen, the description of the Q = 2-quasiparticle by perturbation theory 
breaks down as the chiral angles approach the superintegrable value = (p — |. This, 

a 



of course, is related to the overlapping of the Q = 2-particle with the continuum. At 
<j> = ip = |-, where for P = the Q = 2-particle sits just at the threshold of the two 
Q = 1-particle scattering states (sometimes expressed in the form that for these chiral 
angles m<i = 2ra\ ||). Since one can guess the correct diagonal linear combination of 



states at this point for < A < 1, degenerate perturbation theory is possible here |13 
However, it is not known how to generalize this at = tp = | to P ^ 0, because in the 
standard basis the diagonal combination of states probably is highly complicated. 

Finite-size methods are not affected by these problems and still give good information 
(these will be useful even for (p > |, > |) as we shall see by comparison to the 
exact N — > oo formula (p~2|) . In particular, we are going to check now whether we can 
confirm the details of the peculiar energy-momentum rule of the Q = 2-quasiparticle 
(|P^) by finite-size calculations. After showing that this is possible, we go ahead and look 
whether a similar unusual energy-momentum rule is also found for non-superintegrable 
chiral angles. The question whether integrability plays a special role here will be studied 
choosing chiral angles for which the system probably is not integrable. 

Fig. 4 shows the curves calculated from (|T2| ) together with the 9 lowest N = 12 
finite-size eigenvalues for each Q, P for <p = V 9 = f an d A = 0.50 which is right in the 
center of the high-temperature phase. 

In the lower part of Fig. 4 we see the curve e r (P) drawn according to the analytic 
formula (0). The Q — 1, N = 12-data (crosses) lie within < 2 ■ 10~ 5 on this curve Q. 

The small dashed lines mark the two-particle band expected for two Q = 1-quasipar- 
ticles according to both flT5| ) (rounded lower curve at P ~ 1.37r) and the approximation 
(|T7|) (crossing curves). For P < 3n/2 the Q = 2-particle dips into the two-particle band, 
which for finite- N (consider only the Q = 2-squares in the Figure) contains many states. 
Studying their iV-dependence one finds that in general these points move considerably 
when N is varied, indicating an only power-convergence to the thermodynamic limit. 
However, one notices that e<i s neatly passes through several N = 12-finite-size eigenvalues 
which for P < Att/3 are some of the higher lying levels. We are now going to argue that it 
is no random coincidence that the e2 S — curve passes through the iV = 12-points but this 
comes because these states belong to an exponentially converging sequence, characteristic 
of a fundamental quasiparticle (For the concept of a quasiparticle see the beginning of 
Sec. 2). 

Exponential convergence for states of fixed momentum P can be inferred from finite- 
size data directly only if a sequence of different iV-values at the same P is accessible. 
Since P = 2nk/N, with k integer, this is possible e.g. for P = ir, where we can use data 
for all even N. 

However, being limited to iV < 12, there are only a few other values of P where 
more than one iV can be used, e.g. P = 2tt/3 which is reached for iV = 12, 9, 6, 3 with 
k = 4, 3, 2, 1. Table 1 gives such data for the case A = 0.5 of of Fig. 4. We compare 
the convergence in N of the high-lying Q = 2-levels which are very close to the curve 
C2s with the lowest Q = 1-data. About the Q = 1-levels we are on quite safe grounds 
to assume that these approximate the well-isolated Q = 1-quasiparticle exponentially 
fast for N — > oo. Table 1 shows clearly that the convergence in N to the high- lying 
Q = 2-levels is practically as good as that of the Q = 1-levels and consequently also 



2 Also the Ei(P)-third order perturbation curve (0) (not shown in the Figure) agrees still within < 2% 
with the exact expression. Observe that the A here is much smaller than the A = 0.9 chosen in Fig. 1. 
Compare also the fourth columns in Tables 1 and 2 below. 
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Table 1: Comparison of the exact formulae (T2) with finite chain size data for iV = 12 
and less sites and with the 3ro?-order perturbation formula for A = 0.50 and <fi = ip = | 
(superintegrable case). In the last column we indicate how high up in its (Q = 2,P)- 
sector the quoted energy level appears. For Q — 1 the finite-iV-values are always the 
Q = 1-ground state (1st) levels. 



exponential. This way we can trace the Q = 2-quasiparticle, if it is present, over its 
whole range of the momentum by our finite-size methods. Although at other values of 
P (e.g. P = 7tt/6) we have only one "measurement" for N = 12, we conclude from 
the closeness of the N = 12-value to the exact result that also here we probably have 
exponential convergence. In contrast, in the case of power convergence, typically the 
N = 12- values are off the N — > oo-limit still by several percent. 

We like to emphasize that the choice of levels, even if we get very high up in the spec- 
trum, in practically all cases is quite unambiguous: e.g. for A = 0.5 and P = -y the next 
neighbouring levels to those quoted in Table 1 are: for N = 12: 5.8192/6.0000/6.1965, 
for N = 9: 5.6962/5.9998/6.2942 and for N = Q: 5.3865/5.9964/7.8306. For A = 0.5 
and P = 7r at iV = 12 we have the neighbours 5.2154/5.4755/5.9303 which are also well 
separated. That we are not victims of coincidences will be very clear from Fig. 5 below 
and also from Table 2 which gives the analogous results for A = 0.90. Here, since we 
are much closer to the phase transition line, the level structure is much more involved. 
Nevertheless, although the mass scale of the system is smaller by a factor 5 as compared 
to A = 0.5 (as a relevant scale we may consider Ei(P = 0) = 2(1 — A)), and we have to 
go up to the 26th level at P = 2tt/3 and iV = 12, the convergence is still excellent. 

Notice that the convergence is not monotonous in N but with oscillations superim- 
posed, a phenomenon quite common in the chiral Potts model ||. 

To show that below P = 2ir/3 there is no exponentially converging level is not so easy 
because we do not know precisely where to search. We cannot search in the immediate 
neighbourhood below P = 2n/3, because for N < 13 sites (with some effort, iV = 13- 
chains can be calculated but are quite useless here since 13 is a prime number), the next 
lowest point where we can test for fast convergence is P = ir/2. There we do not see 
any convincingly close levels in the N = 12, k = 3 and N = 8, k = 2-spectra. A better 
argument that P = 2it/3 is the lower endpoint of the Q = 2-quasiparticle with respect 
to P will be given below in connection with Fig. 4, using the change of the spectrum at 
P = 2tt/3 with 0, ip when moving away from the superintegrable line. 

Summarizing, we see that even from our numerical studies one obtains strong evidence 
that one can interpret the spectrum of ([!]) on the superintegrable line in terms of two 
elementary quasiparticles with energies e r (P) and e2 S {P)- The unusual features of the 
Q = 2 quasiparticle are clearly reflected in the finite-size data. Decreasing P from 2tt 
towards lower values, the Q = 2 excitation crosses the scattering threshold of two Q = 1 
excitations at some point but stays stable. 



4.2 The non-superintegrable case 

We expect that the limitation of the Q = 2-quasiparticle to only part of the range of 
the momentum should not only appear at the superintegrable line (f> — ip = | but by 
continuity should be valid, maybe with changing range, also for neighbouring values of 
(f) and if. 

From our earlier calculations in |TD| we get a hint about the endpoint at P = 2n\ 



For = if > -| this endpoint should move to P < 2n because we found that the 
Q = 2— particle becomes unstable at P = (or P = 2ti) (which was the only value for 
P considered there). 

In order to get more information, we now study the Q = 2 spectra at P = 2ir/3 for 
iV = 6, 9, 12 sites as functions of the chiral angle <fi at fixed A = 0.50. In order to vary 
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Table 2: Comparison of the exact formulae ( flj3 ) with finite chain size data for iV < 12 
sites and with the 3r<i-order perturbation formula for A = 0.90 and = tp = | (super- 
integrable case). In the last column we indicate how high up in its (Q = 2,P)-sector 
the particular quoted energy level appears. In contrast to Table 1, here for Q = 1 the 
levels approximating the curve Ei(P) are not always the Q = 1-ground state levels. 
Recall Fig. 1 which showed that the Q = 1— quasiparticle is entering the continuum for 

27T/3 < P < 7T. 
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only one parameter, we have used both the self-dual choice = ip and the integrable 
choice (H) but we found that both behave very similarly. So here we present only results 
for the self-dual case. For this choice ([!]) probably is not integrable, and the effects we 
are going to see should not be due to the integrability. 

In Fig. 5 one clearly distinguishes three curves, one for each N = 6, 9, 12, which 
follow each other very closely for > 90° and separate in the region < 90°. At = 90° 
this is the Q = 2— particle of (|T2"D, and the three values for N = 6, 9, 12 are those of 
Table 1 (5.9999978 etc.). Notice that indeed there are no other levels close to 6.000 
around = 90°, and that nowhere else in the picture are lines from N = 6, 9, 12 
following closely the same course. We see that the Q = 2-quasiparticle curve can be 
followed uniquely up to = 125°, and that the fast convergence from N = 6 to N = 12 
is still there. Around ~ 120° the convergence is quite perfect. The energy values at 
= 120° being 5.57251; 5.57211; 5.57192 for iV = 6, 9, 12, respectively. 

For larger angles the Q = 2-quasiparticle level gets so high up in the spectrum that 
computational problems limit to follow it for N = 12: at = 126° it is already the 30th 
level in its sector. In Fig. 5 we have plotted up to the 29th N = 12-level in the Q = 2, 
P = 27r/3-sector. Around ~ 100° the convergence is somewhat worse but still quite 
good: 5.81200; 5.79749; 5.79150 at = 99°. We interpret this as continuing presence 
and stability of the Q = 2-quasiparticle for 90° < < 125° at P = 2vr/3. 

However, for < 90° the curves for N = Q, 9, 12 clearly separate. We consider this 
to be an indication that here the Q = 2— particle gets out of its stability range. This is 
consistent with the analytic result that for = (p — |, P — 2-7r/3 is the lower momentum 
bound of the second quasiparticle. 

On the scale on which Fig. 5 is drawn it is difficult to show that apparent crossings 
of level-lines with the same N, Q and P are in fact avoided crossings, since the levels 
come very close to each other. E.g. consider the N = 9-levels which around ~ 99.035° 
seem to cross at E « 5.797. At their closest encounter they keep separated in energy by 
AE « 0.0007911. 

It remains to see what happens above the superintegrable line for other momenta. 
E.g. for P = 7r we find that the three Q = 2-levels for N = 8, 10, 12, which were very 
close to each other for = <p = |, stay very close together if we move in 0. Table 3 gives 
a few numbers for the energy of the Q = 2-particle obtained by following this triple. 
At intermediate momentum values, where we have only one N < 12 available, we can 
make use of a property which is clearly seen in Fig. 5 but which is qualitatively the same 
for other P, too: the levels approximating the Q = 2-quasiparticle have a slope d^/difi 
which is very different from that of other levels. Using this feature we can follow the 
quasiparticle line to larger 0, starting from the known value at = tp = |. Of course, 
when using just this method we have no direct control over the size of errors. 

In order to illustrate that the integrable case is analogous to the self-dual one for the 
parameters under consideration we also mention two further values of ^{P): For the 
integrable choice of parameters A = 0.50, <p = 114° one finds E2(P = 27r/3) = 5.683(2) 
and E.2{P = tt) = 5.30608(3). Note the similarity between these two values including 
their errors with the corresponding ones in Table 3. 

We conclude this Section with the remark that for small A the behaviours of the 
self-dual (0 = <p) and the integrable versions (|4]) have to be different, because in the 
integrable case for small A the non-hermitian region comes close to the line = ip = |, 
whereas in the self-dual case little special is happening for small A when moving in 
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= <P 


90° 


96° 102° 


108° 


114° 


120° 


P = 2tt/3 

P = 71 


6.0000(1) 
5.475(1) 


5.84(1) 5.745(5) 
5.4261(5) 5.3827(2) 


5.687(7) 
5.34024(6) 


5.626(2) 
5.29667(4) 


5.5718(3) 
5.25038(5) 



Table 3: Q = 2-quasiparticle energies at A = 0.50 for two values of the momentum P 
and varying chiral angles = (p. The values in the Table are estimates for the limit 
iV — > oo obtained from the values at iV = 12, 9, 6 and N = 12, 10, 8 sites, respectively. 



around <j> = ~ . 



5 Correlation functions 

In this Section we define and study a non-translationally invariant correlation function 
for the operator V of the Z 3 -chain numerically. We shall use the definition 

Cr{x) := <" ' fo r ' ' »> < x < £ (19) 

{v I v) 2 

where | v) is the lowest eigenvector of the Hamiltonian in the sector Q = P = which 



in the massive high-temperature phase is the ground state. The definition (0^) ensures 
Cp(— x) = Cr{x)*. Non-translationally invariant correlation functions are difficult to 
measure experimentally but for our purpose it will be crucial to work really with r+ +1 ri 
and not to replace it by N' 1 Y.r=i ^t+r^ri because this would destroy the oscillation. 

Because we are studying a massive phase, we expect an exponential decay in addition 
to the oscillatory contribution. So a simple ansatz for Cr{x) may be made in terms of a 
complex correlation length or equivalently, in terms of the real correlation length £r and 
the oscillation length L defined through 

C T {x) = ae 2mx/L e- x '^ + (1 - a)8 Xj0 . (20) 

For the massive low-temperature phase a non-vanishing wave vector has been pre- 
dicted in H where also its critical exponent was calculated from level crossings. Per- 
turbative calculations of the correlation function Cr(x) for the massive low-temperature 
phase have been performed in |TT| and indications for an oscillatory contribution have 



been obtained. Detailed perturbative studies of Cr(x) in the massive high-temperature 
phase will be presented in fll3| . 

We have noticed in |TU[ that in the high-temperature limit, the dispersion curve of 



the Q = 1-particle has its minimum at P m i n = Re0/3, as it is obvious from the term of 
first order in A of Ei(P), see eq. (0). The value of P m i n decreases with increasing A: e.g. 
for (f) = ip = tt/2 one finds from (0), QTJ) that at A = 0.901292 P min w 0.225 Q so that 



at this A, P m in has diminished by a factor of w 2.33 from its value at small A. In Sec. 2 
(see Fig. 1) we have also noticed that for large A the third order formula (0) does not 
describe the dispersion curve around the minimum. 

3 More precise parameter values for the transition to the IC phase for <f) = cp = ^ are: 
A = 0.90129284..., P min = 0.22503375(3). 
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The non-zero shift of the minimum of the dispersion curve has the result that at 
7^ 0, already for A — > 0, the first excited state is no longer at P = 0. For a finite 
iV-site chain with chain momentum k = NP/2n, the lowest gap has Ak = [-/V\P min /27r]. 
This gives rise to a proliferation of cross-overs in the first excited states with increasing 



iV which in turn, via finite-size scaling arguments [15], indicates an oscillating AP ^ 
correlation function. So we may expect P m i n to be related to the oscillation length L of 
the correlation function. 

For A —>■ this relation is readily found using a high-temperature expansion for the 
correlation function flHf). Calculating ( [19] ) up to second order in A, we obtain 

i0/3 \2 

Cr(x) = S x , + XS xl j- + - e~ 2 ^ 3 + 2S X , 2 e 2 ^ 3 } + 0(A 3 ). (21) 

3 cos (|J 12 cos (|j L J 

In order to extract the complex correlation length from this formula to lowest nontrivial 
order, we consider 

C r (2) e^/ 3 



C r (l) 2cos(f)" 
On the other hand, from (|20|) we have 



(22) 



Cr(l) 

so that, comparing (P2"|) with (123) we obtain 



Cr(2) -expl^-f I, (23) 




L Re0 = 6tt, f£ 1 = ^ - In — * (24) 

3 2cos(i|) 



(We write Re0, since e.g. in the integrable case (|J) can be complex for small A). Using 
the small- A expression for P m m, the first of these equations can be rewritten as a relation 
between the oscillation length L and the minimum P m i n of the dispersion relation (|^): 

P min L = 2vr, (25) 

which is valid at general 0, if and small inverse temperatures A. 

In the following, we will show that with a finite-size numerical calculation, using up 
to 13 sites, one can get useful information on the oscillation length L at various values of 
the parameters in the high-temperature phase. Since, as we have seen, L is conceptually 
linked to P m i n , we shall conveniently give our results in terms of the product P m i n L. It 
will turn out that ( plf) seems to be satisfied approximately for a wide range of parameters. 

In order to evaluate the correlation function ( |l9|) numerically, one needs the ground 
state of the Hamiltonian ([I]). The ground state of the Z3-chain is easily obtained for 
general values of the parameters (even in the non-hermitian case G (C [14[]) using vector 
iteration up to iV = 13 sites. The more difficult point is to calculate the matrix elements 
of r^ +i r 1 , because this operator does not conserve momentum and thus does not leave a 
space of momentum eigenstates invariant. Table 4 shows the correlation function C-p(x) 
obtained in this manner for iV = 12 and iV = 13 sites at two different points in the 
phase diagram: first for A = 0.5 in the superintegrable case, and second, for A = 0.25 
and the chiral angle if above the superintegrable line and chosen complex (so that the 
Hamiltonian becomes non-hermitian) satisfying the integrability condition @. 



1 o 







A = 0.50 


<p=&, = tt + 1.70005 i, A = 0.25 
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.00028 - .00052i 


6 
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-.00063 + .00058? 


.00037 


.00022 - .00005i 



Table 4: Numerical results for the correlation function C?{x) at iV = 12, 13 sites 



The excellent agreement between the calculations at iV = 12 and iV = 13 sites shows 
that the contribution of boundary terms to the correlation functions is small. So we can 
use the N = 12— values as a good approximation for the infinite chain limit as long as 
x < 6. 

We thus obtain six complex numerical values for Cr at each fixed point (A, <p, <p) in 
the phase diagram considered. By a fitting procedure, we want to translate these six 
values into the three parameters £r> L and a of (pop. We proceed as follows: First, we 
obtain £ r by calculating Re(ln(Cr(x)/Cr(x + l))) -1 and averaging over x. Next, the first 
zero of Ke(e x ^ T Cr(x)) is estimated by linear interpolation for two neighbouring values 
and L/4 is obtained by averaging. Finally, a is chosen such that the difference 

ReC r (x) - ae~ xliv cos (2ttx/L) (26) 

is minimal for x = 1,2. Table 5 collects our results for 12 choices of A and the chiral 
angles. 

For the entries marked with a '*' ReCr does not yet reach zero in the interval [0,6]. 
Here we calculated L from the formula 

t = I v — m\ 

4^ 1 Im(ln(C r (x)/C r (x + l)))- 1 1 

The values of P m i n in Table 5 have been obtained by first calculating Ei(P) numerically 
at iV = 12 sites and afterwards minimizing the finite Fourier decomposition of (Ei(P)) 2 
numerically - compare also Table 8 of JjlJ. These values are approximations with errors 
which are difficult to estimate. 

The eighth and ninth line in Table 5 give situations close to, and in the IC phase. In 
both cases Ei(P) partly enters two-particle scattering states. For = ^ = ^, A = |we 
are in the IC-phase with Ei(P) becoming negative. We nevertheless define the correlation 
function for these cases with | v) still denoting the Q = P = 0— ground state, although it 
is no longer the lowest state of (|l|). Since in these cases the method explained above for 
estimating P m i n becomes impractical, here we proceeded differently: We first determine 



1 





A 




L 


a 


p min 


-^min / 2 7T 


3tt 
8 


3vr 
8 


i 

4 


0.53(3) 


18(2) 


0.65(3) 


0.341 


1.0(1)) 


3tt 
8 


3vr 
8 


1 
2 


0.8(2) 


25(1)* 




0.283 


1.14(6) 


3tt 
8 


3tt 
8 


3 
4 


1.8(6) 


50 ± 23* 




0.211 


1.7(8) 


7T 

2 


7T 

2 


1 
4 


0.6(1) 


14.05(6) 


0.44(8) 


0.470 


1.051(4) 


7T 

2 


7T 

2 


1 
2 


0.60(5) 


16(1) 


0.59(3) 


0.401 


1.07(8) 


7T 

2 


7T 

2 


3 
4 


1.3(2) 


29(6)* 




0.189 


0.9(2) 


3tt 
4 


3vr 
4 


1 
4 


0.62(8) 


10(1) 


0.4(2) 


0.746 


1.18 ±0.12 


3tt 
4 


3tt 
4 


1 
2 


1.2(6) 


12(1) 


0.4(1) 


o.esst 


1.34 ±0.14 


3tt 
4 


3vr 
4 


3 
4 


1.4(1) 


16(1) 


0.65(3) 


0.602* 


1.54(8) 


3tt 
8 


-0.98942i 


1 
4 


0.7(2) 


oo 









3tt 
8 


0.69919 


1 
2 


1.1(3) 


49(7)* 




0.159 


1.2(2) 


3tt 
4 


7r + 1.70004? 


1 
4 


1.0(5) 


5.9(1) 


0.3(1) 


7T 

3 


0.99(2) 



Table 5: Parameters for the correlation function (|20"D calculated numerically for iV = 12 
sites. In the first nine lines we give data for the choice (j) — tp of the chiral angles, where 
three lines correspond to the superintegrable case <p — V 9 = f • The bottom three lines 
use the integrable choice (|j). 

the smallest energy gap in the Q = 1 sector at iV = 8, . . . , 12 which appears at the 
smallest value possible for the discrete momentum, P = 2n/N. Next, a polynomial 
interpolation between these five values for the energy gap was minimized numerically. 
In Table 5 we marked these results by 

At the end of Table 5 we give three examples with integrable choices of the chiral 
angles (f|). For ip = A = |, the integrable Hamiltonian is still hermitian and behaves 
very similarly to the the self-dual case. However, for small values of A and ^ 7r/2, 
one has to make one of the angles (p, <ft complex in order to satisfy (f|). We choose <fi 
complex. Below the superintegrable line <p = p = | one can take (p purely imaginary, 
whereas above the superintegrable line one has to admit a non-vanishing real part, e.g. 
(j) — Ti G i R . With these choices one can verify, using e.g. (0), that Ei(— P) = Ej.(P)* for 
(p < | and imaginary (f> (for details, see |TJJ). Therefore, it seems sensible to set P m i n = 0, 
and, if ( p5|) holds at least approximately, we expect no oscillation. Indeed, the complete 
correlation function C*r(x) is real for the integrable choice of parameters ip = A = | 
which in view of (PC|), implies L = oo. For the integrable point <p — ^f, A = \ above 
the superintegrable line we have verified that (|7|) and a numerical evaluation of Ei(P) 
at JV = 12 sites satisfy Ei(-P - | ) = Ei(P - |)*. Thus, we set P min = f which gives 
excellent agreement with the prediction LP min = 2it. 

That the above procedures yield reasonable fits is demonstrated by Fig. 6 which 
shows the stretched correlation function e x ^ r Cr(x) obtained numerically for iV = 13 
sites at the integrable point p = ^f, A = -r in comparison to the fits. The errors have 
been estimated from the difference between the results for iV = 12 and iV = 13 sites in 
Table 4. The agreement for all x, not only in the real part but also in the imaginary 
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part is good. In particular, the oscillation is clearly visible. 

The conjecture LP m i n = 2tt is well verified for all values in Table 5, bearing in mind 
that we have ignored systematic errors. This may apply also to the point = <p = ^f, A = 
| in the massless IC phase, if we had underestimated the error in the determination of 
L. 



6 Conclusions 

In this paper we have presented further numerical evidence that the spectrum of the 
general Z 3 -spin quantum chain (|I|) in the massive high-temperature phase can be in- 
terpreted in terms of two fundamental quasiparticles and their scattering states. Using 
duality 0] our results about the quasiparticle spectra can be pulled over to the massive 



low-temperature phase. 

The quasiparticle carrying Z 3 -charge Q = 2 shows a quite peculiar energy- mo men- 
turn relation. It becomes unstable for a certain range of the momentum, while exhibit- 
ing a remarkable stability against decay into multi-particle states allowed by energy- 
momentum conservation in other ranges of the momentum. Since we find this feature 
not only in the superintegrable case, this generalizes the analytic result found by the 
Stony-Brook group || to cases for which no exact methods are known. 

We also calculated the correlation function for the operator T using a numerical 
evaluation of the ground state of the Z 3 -model for up to 13 sites. Although this approach 
is limited to short ranges, we are able to estimate correlation lengths in the massive high- 
temperature phase and can show that the correlation functions oscillate. We find that 
the product of the oscillation length L and the minimum of the dispersion curve P m i n 
(which is non-zero due to parity violation) satisfies LP m \ n « 27r. 
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Figure captions 



Fig. 1: Momentum dependence of the spectrum of the superintegrable Z3-Hamiltonian 
in the Q = 1-sector for A = 0.90. Symbols: Finite chain data, dashed curve: exact 
lowest level from 0, point-dashed curve: 3rd-order perturbation in A according to 
eq. (0). 

Fig. 2: Spectrum of the Z 3 -Hamiltonian at <f> = ip — ?, A = 0.10. All three charge 
sectors have been plotted in the same figure. The dots indicate eigenvalues for iV = 
12 sites whereas the lines are obtained from the 3r(i-order perturbative expansions 
(0) and (|10D . We have limited ourselves to show for each (Q,P) the nine lowest 
eigenvalues only. One observes that the eigenvalues nicely group within the band 
boundaries expected according to fllSf) . For Q — 1 and Q = 2 the first levels 
building up the three-particle bands are seen for E > 8.5. For a discussion of the 
four Q = 0-points at P = n outside the (Ei + E 2 )-band, see the text. 

Fig. 3: Same as Fig. 2 but for (p = p = A = 0.25. Note that all three charge sectors 
have been plotted in the same figure. Thus, seemingly overlapping bands belong 
to different charge sectors and so are no real level crossings. The dots indicate 
eigenvalues for N = 4, . . . , 12 sites. The curves show the perturbative results for 
Ei(P) (full curve), Eq =2 (P) (see eq. ([15]), long dashes) and E 2 (P) (fine dashes), 
according to (|7|), ([TP]) . 

For some range of P the Q = 2 excitation enters the energy band of two Q = 1 
particles. In this region, the Q = 2 "particle" is difficult to trace, although outside 
this region the convergence in N is very fast. This indicates that the Q = 2 
fundamental excitation is not stable for all momenta P. 

Fig. 4: Spectrum of the superintegrable Z 3 -Hamiltonian &t<j) = p = ^,\ = 0.5. All 
three charge sectors have been plotted in the same figure but only iV = 12-finite- 
size data are shown. The dashed curves show the analytic result of [0] , eq. ( fT2] ) for 
2Ei(P/2) (long dashes) and Eq =2 (P) (fine dashes). The full curve is the analytic 
result for E 2 (P). Note that the Q = 2 fundamental particle exists only for P > 

Fig. 5: Part of the spectrum in the Q = 2-sector at A = 0.50 and fixed momentum 
P = 2tt/3 for N = 6,9, 12 sites, as a function of the chiral angle <p = <p (self-dual 
choice). On sees that just one level from iV = 6,9,12 each moves very closely 
together with their companions from the other iV when <p varies. We take this 
as evidence for an exponential convergence of this level in N , characteristic for a 
single-particle state. In the upper right corner there are more iV = 12-levels which 
have not been drawn since we limited ourselves to calculate for iV = 12 no more 
than 29 levels. 

Fig. 6: Numerical values for the correlation function Cr{x) at iV = 13 sites stretched 
by e x/ tr in comparison to the fit (^) at = ^f, A = \, <p = ir + 1.700041 The 
error bars are given by the differences between N = 12 and iV = 13 sites. The 
parameters used for the fit are £r = 0.77, L = 6 and a = 0.35. One observes that 
the agreement with these fits is good and the oscillatory contribution to Cr(x) is 
clearly visible. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9409170vl 



Fig. 1: Spectrum of the Q=l sector at (]) = (p = 7C/2, X = 9/10 
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This figure "figl-2.png" is available in "png" format from: 
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Fig. 2: Spectrum at § = (p = rc/4, X = 1/10 
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This figure "figl-3.png" is available in "png" format from: 
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Fig. 3: Spectrum at § = (p = 3 7i/8, X = 1/4 
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Fig. 5: Spectrum of the Q=2 sector for X = 1/2, P=2 7t/3 
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Fig. 6: Correlation for T at the integrable point (p = 3 7C/4, A, = 1/4, (|) = 71 + i 1.70004 



0.4 



0.2 



^ 



-0.2 



-0.4 



-0.6 







X 



X 



exp(x/^) Re(C r (x)) h> 
exp(x/^) Im(C r (x)) 
0.35 cos(2 71 x/6) 
0.35 sin(2 71 x/6) 



h-x — i 



<y 



x 



<> 



